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The case is strong that cosmic microwave background (CMB) and large scale 
structure (LSS) observations can be combined to determine the theory of struc- 
ture formation and the cosmological parameters that define it. We review: the 
relevant 10+ parameters associated with the inflation model of fluctuation gener- 
ation and the matter content of the Universe; the relation between LSS and pri- 
mary and secondary CMB anisotropy probes as a function of wavenumber; how 
COBE constraints on energy injection rule out explosions as a dominant source 
of LSS; and how current anisotropy band-powers in multipole-space, at levels 
~ (10 -5 ) 2 , strongly support the gravitational instability theory and suggest the 
universe could not have reionized too early. We use Bayesian analysis methods 
to determine what current CMB and CMB+LSS data imply for inflation-based 
Gaussian fluctuations in tilted ACDM, AhCDM and oCDM model sequences with 
cosmological age 11-15 Gyr, consisting of mixtures of baryons, cold "c" (and pos- 
sibly hot "h" ) dark matter, vacuum energy "A" , and curvature energy "o" in open 
cosmologies. For example, we find the slope of the initial spectrum is within about 
5% of the (preferred) scale invariant form when just the CMB data is used, and 
for ACDM when LSS data is combined with CMB; with both, a nonzero value of 
is strongly preferred (~ 2/3 for a 13 Gyr sequence, similar to the value from 
SNIa). The oCDM sequence prefers Qtot < 1 ; but is overall much less likely than 
the flat Q\ ^ sequence with CMB+LSS. We also review the rosy forecasts of 
angular power spectra and parameter estimates from future balloon and satellite 
experiments when foreground and systematic effects are ignored to show where 
cosmic parameter determination can go with just CMB information alone. 



1. The Relation Between CMB and LSS Observables 

In this section, we first present an overview of the relation between the scales 
that CMB anisotropies probe, those that large scale structure observations of 
galaxy clustering probe, and the scales that are responsible for collapsed object 
formation in hierarchical models of structure formation, in particular those deter- 
mining the abundances of clusters and galaxies. We review the basic parameters 
of amplitude and tilt characterizing the fluctuations in the simplest versions of 
inflation, but consider progressively more baroque inflation models needing pro- 
gressively more functional freedom in describing post-inflation fluctuation spec- 
tra. We then describe the high precision that has been achieved in calculations of 
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primary CMB anisotropies (those determinable with linear perturbation theory), 
and the less precisely calculable secondary anisotropies arising from nonlinear 
processes in the medium. 

(a) CMB as a Probe of Early Universe Physics 

The source of fluctuations to input into the cosmic structure formation problem 
is likely to be found in early universe physics. We want to measure the CMB (and 
large scale structure) response to these initial fluctuations. The goal is to peer into 
the physical mechanism by which the fluctuations were generated. The contenders 
for generation mechanism are (1) "zero point" quantum noise in scalar and tensor 
fields that must be there in the early universe if quantum mechanics is applicable 
and (2) topological defects which may arise in the inevitable phase transitions 
expected in the early universe. 

From CMB and LSS observations we hope to learn: the statistics of the fluc- 
tuations, whether Gaussian or non-Gaussian; the mode, whether adiabatic or 
isocurvature scalar perturbations, and whether there is a significant component 
in gravitational wave tensor perturbations; the power spectra for these modes, 
V$(k),Vi s (k),VGw(k) as a function of comoving wavenumber k = 2ira/\, with 
the cosmological scale factor a(t) removed from the physical wavelength X(t) and 
set to unity now. The length unit is h _1 Mpc, where h is the Hubble parameter 
in units of 100 km s _1 Mpc -1 , i.e., really a velocity unit. Until a few years ago 
h was considered to be uncertain by a factor of two or so, but is now thought to 
be between 0.6 and 0.7. 

Sample initial and evolved power spectra for the gravitational potential V$(k) 
(= da\/d\n k, the rms power per dink band) are shown in Fig. |l[ The (linear) 
density power spectra, V p {k) oc k A V$(k), are also shown in Fig. [j]. (We use 
V(k) = k s P(k)/(2ir 2 ) for power spectra, the variance in the fluctuation variable 
per Ink, rather than the oft-plotted mean-squared fluctuation for mode k, P(k), 
so V p = A 2 (k) in the notation of Peacock (1997).) As the Universe evolves the 
initial shape of V§ (nearly flat or scale invariant) is modified by characteristic 
scales imprinted on it that reflect the values of cosmological parameters such 
as the energy densities of baryons, cold and hot dark matter, in the vacuum 
(cosmological constant), and in curvature. Many observables can be expressed as 
weighted integrals over k of the power spectra and thus can probe both density 
parameters and initial fluctuation parameters. 

(b) Cosmic Structure and the Nonlinear Wavenumber 
In hierarchical structure formation models such as those considered here, as the 

1/2 

universe evolves Vp (k) grows with time until it crosses unity at small scales, and 
the first star forming tiny dwarf galaxies appear ("1st *"), typically at a redshift 
of about 20. The nonlinear wavenumber fc/vx(£) ; defined by Jq NL V p (k)d\nk=\, 
decreases as the universe expands, leaving in its wake dwarf galaxies (dG), normal 
galaxies (gal), groups (gps) and clusters (els), forming from waves concentrated 
in the A:-space bands that their labels cover in Fig. [I]. Equivalent mass scales are 
given above them. 

Scales just below kjsn are weakly nonlinear and define the characteristic pat- 
terns of filaments connecting clusters, and membranes connecting filaments. Voids 
are rare density minima which have opened up by gravitational dynamics and 
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merged, opposite to the equally abundant rare density maxima, the clusters, in 
which the space collapses by factors of 5-10 and more. 

At k > k]y i (t) , nonlinearities and complications associated with dissipative gas 
processes can obscure the direct connection to the early universe physics. Most 
easily interpretable are observables probing the linear regime now, k < fc/vx(io)- 
CMB anisotropies arising from the linear regime are termed primary. As Fig. [l] 
shows, these probe three decades in wavenumber, with the high k cutoff defined 
by the physics at z ~ 1000 when CMB photons decoupled, not k^L at that time. 
Within the LSS band, two important scales for the CMB arise: the sound crossing 
distance at photon decoupling, ~ 100 h _1 Mpc, and the width of the region over 
which this decoupling occurs, which is about a factor of 10 smaller, and below 
which the primary CMB anisotropies are damped. LSS observations of galaxy 
clustering at low redshift probe a smaller range, but which overlaps the CMB 
range. We have hope that z ~ 3 LSS observations, when fc/vz,(£) was larger, can 
extend the range, but gas dynamics can modify the relation between observable 
and power spectrum in complex ways. Although probes based on catalogues of 
high redshift galaxies and quasars, and on quasar absorption lines from the in- 
tergalactic medium, represent a very exciting observational frontier, it will be 
difficult for theoretical conclusions about the early universe and the underlying 
fluctuations to be divorced from these "gastrophysical" complications. Secondary 
anisotropies of the CMB (§ |h~|), those associated with nonlinear phenomena, also 
probe smaller scales and the "gastrophysical" realm. 

(c) Probing Wavenumber Bands with the CMB and LSS 

Although the scales we can probe most effectively are smaller than the size of 
our Hubble patch (~ 3000 h _1 Mpc), because ultralong waves contribute gentle 
gradients to CMB observables, we can in fact place useful constraints on the 
ultralarge scale structure (ULSS) realm "beyond our horizon". Indeed current 
constraints on the size of the universe arise partly from this region and partly 
from the very large scale structure (VLSS) region. (For compact spatial manifolds, 
the wavenumbers have an initially discrete spectrum, and are missing ultralong 
waves, limited by the size of the manifold.) 

The COBE data and CMB experiments with somewhat higher resolution probe 
the VLSS region very well. Density fluctuations are highly linear in that regime, 
which is what makes it so simple to analyze. One of the most interesting realms 
is the LSS one, in which CMB observations probe exactly the scales that LSS 
redshift surveys probe. The density fluctuations are linear to weakly nonlinear 
in this realm, so we can still interpret the LSS observations reasonably well - 
with one important caveat: Galaxies form and shine through complex nonlinear 
dissipative processes, so how they are distributed may be rather different than 
how the total mass is distributed. The evidence so far is consistent with this 
"bias" being only a linear amplifier of the mass fluctuations on large scales, albeit 
a different one for different galaxy types. Detailed comparison of the very large 
CMB and LSS redshift survey results we will get over the next five years should 
help enormously in determining the statistical nature of the bias. 

Because the V§ of the COBE- normalized sCDM model shown shoots high rel- 
ative to the cluster data point, the sCDM model is strongly ruled out. More 
rigorous discussion of what is compatible with COBE, smaller angle CMB exper- 
iments such as SK95, the cluster data point and the shape of the V$ spectrum as 
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estimated from galaxy clustering data is given in § |dj The filter functions plotted 
for SK95, Planck, etc. show the bands they are sensitive to: multiplying by a 
fe-space AT/T power spectrum gives the variance per In A; (e.g., Bond 1996). 

(d) The Cosmic Parameters of Structure Formation Theories 

Even simple Gaussian inflation-generated fluctuations for structure formation 
have a large number of early universe parameters we would wish to determine 
(§ E]) : power spectrum amplitudes at some normalization wavenumber k n for 
the modes present, {7 3 $(&n))'Pis(&n)>'pGW'(&n)}; shape functions for the "tilts" 
{u s (k),Ui s (k),Ut(k)}, usually chosen to be constant or with a logarithmic cor- 
rection, e.g., v s {k n ), dv s (k n )/d\n k. (The scalar tilt for adiabatic fluctuations, 
v s (k) = dlnV^/dlnk, is related to the usual index, n s , by v s = n s — 1.) The 
transport problem (§ |g]) is dependent upon physical processes, and hence on 
physical parameters. A partial list includes the Hubble parameter h, various 
mean energy densities {&tot, ^cdm, ^/i<2m}h 2 , and parameters character- 

izing the ionization history of the Universe, e.g., the Compton optical depth tc 
from a reheating redshift z re h to the present. Instead of Qtoti we prefer to use the 
curvature energy parameter, = 1 — Qtot, thus zero for the flat case. In this 
space, the Hubble parameter, h = (^(^jh 2 )) 1 / 2 , an d the age of the Universe, 
to, are functions of the fijh 2 . The density in nonrelativistic (clustering) particles 
is Q nr = + ^cdm + ^hdm- The density in relativistic particles, Q er , includes 
photons, relativistic neutrinos and decaying particle products, if any. f2 er , the 
abundance of primordial helium, etc. should also be considered as parameters 
to be determined. The count is thus at least 17, and many more if we do not 
restrict the shape of V$(k) through theoretical considerations of what is "likely" 
in inflation models. Estimates of errors on a smaller 9 parameter inflation set for 
the MAP and Planck satellites are given in § ^J. 

The arena in which CMB theory battles observation is the anisotropy power 
spectrum in multipole space, as in Figs. Qj3|, which show how primary C/s vary 
with some of these cosmic parameters. Here Ce = l{t + l)(|(AT/T)£ m | 2 )/(2-7r). 
The C/s are normalized to the 4-year <imr(53+90+31)(A+B) data (Bennett et 
al. 1996a, Bond 1995, Bond h Jaffe 1997). The arena for LSS theory battling 
observations is the V$ of Fig. [l|. (Usually it is V p /k 3 ~ kV§ which is plotted.) 

For a given model, the early universe V<& is uniquely related to late-time power 

1/2 

spectrum measures of relevance for the CMB, such as the quadrupole C 2 or 

1/2 

averages over ^-bands B, (Q)^ , and to LSS measures, such as the rms density 
fluctuation level on the 8h _1 Mpc (cluster) scale, as, so any of these can be used 
in place of the primordial power amplitudes in the parameter set. In inflation, 
the ratio of gravitational wave power to scalar adiabatic power is Vgw/'P® ~ 
— (100/9)^/(1 — ft/2), with small corrections depending upon v s — vt (Bond 
1994, 1996). If such a relationship is assumed, the parameter count is lowered by 
one. 



f It is becoming conventional to refer to f2 nr as Q 
Phil. Trans. R. Soc. Lond. A (1997) 



Constraining LSS with the CMB 5 
0.0001 0.001 0.01 0.1 1 10 100 



1000 
100 
10 



1 



0.1 
0.01 
0.001 



1 

0.5 


0.0001 0.001 0.01 0.1 1 10 100 

comoving wavenumber k hMpc -1 

Figure 1. The bands in comoving wavenumber k probed by CMB primary and secondary 
anisotropy experiments, in particular by the satellites COBE, MAP and Planck, and by large 
scale structure (LSS) observations are contrasted. The width of the CMB photon decoupling 
region and the sound crossing radius (Ar 7 dec, c s r 7( ; ec ) define the effective acoustic peak range 
for primary anistropies (those involving linear fluctuations). Secondary anisotropies arise only 
once matter has gone nonlinear. Sample (linear) gravitational potential power spectra (actually 

Vy 2 {k)) are also plotted, and the y-axis values refer to V^/ 2 /10~ 5 (which is dimensionless) . 
The horizontal dotted line is the post-inflation scale invariant power spectrum, which is bent 
down as the universe evolves by an amount dependent upon the matter content. The hatched 
region at low k gives the 4 year DMR error bar on the $ amplitude in the COBE regime. The 
solid data point in the cluster-band denotes the $ constraint from the abundance of clusters 
(for fit o t = l,f2A=0). The open circles are estimates of the linear $ power from current galaxy 
clustering data by Peacock (1997). A bias is "allowed" to (uniformly) raise the shapes to match 

the observations. The corresponding linear density power spectra, 7-p^ 2 (fc), are also shown rising 
to high k. Models shown in Fig. hf are the "standard" n s = 1 CDM model (labelled T = 0.5 
with Q nr = 1, h = 0.5), a tilted (n s = 0.6, F = 0.5) CDM model and a model with the shape 
modified (r = 0.25) by changing the matter content of the Universe, e.g., Q nr — 0.36, h = 0.7. 
The bands at high k associated with object formation (els, gal, etc. ) and the filters showing the 
bands various CMB experiments probe are discussed in the text. 
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Figure 2. The Ce anisotropy bandpower data for experiments up to summer 1998 are shown in the 
upper left panel. The data are optimally combined into 9 bandpower estimates (with one-sigma 
errors) shown in the upper right panel. To guide the eye an untilted COBE-normalized sCDM 
model is repeated in all panels. The rest of the panels show forecasts of how accurate Ce will 
be determined for this model from balloon and satellite experiments, with parameters given in 
Table 0. 

(e) Fluctuation Freedom in Inflation 

Many variants of the basic inflation theme have been proposed, sometimes with 
radically different consequences for ~ fc 1 -^), and thus for the CMB sky, 

which is used in fact to highly constrain the more baroque models. A rank- 
ordering of inflation possibilities: (1) adiabatic curvature fluctuations with nearly 
uniform scalar tilt over the observable range, slightly more power to large scales 
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Figure 3. The 9 bandpower estimates from current anisotropy data are compared with various 
13 Gyr model sequences: (1) Ho from 50 to 90, SIa, to 0.87, for an untilted ACDM sequence; 
(2) n s from 0.85 to 1.25 for the H = 70 ACDM model (Q A = .66); (3) fish 2 from 0.003 to 
0.05 for the Ho = 70 ACDM model; (4) Ho from 50 to 65, Sl k from to 0.84 for the untilted 
oCDM sequence; (5) the same for a n 3 — 0.9 oCDM sequence, clearly at odds with the data; 
(6) Ho— 50 sequence with neutrino fractions varying from 0.1 to 0.95; (7) shows an isocurvature 
CDM sequence with positive isocurvature tilts ranging from to 0.8; (8) shows that sample 
defect Ce's from Pen, Seljak and Turok (1997) do not fare well compared with the current data; 
Ce's from (Allen et al. 1997) are similar. The bottom right panel is extended to low values to 
show the magnitude of secondary fluctuations from the thermal SZ effect for the ACDM model. 
The kinematic SZ Ce is significantly lower. Dusty emission from early galaxies may lead to high 
signals, but the power is concentrated at higher £, with possibly a weak tail because galaxies 
are correlated extending into the £ < 2000 regime. 



Phil. Trans. R. Soc. Lond. A (1997) 



8 



J.R. Bond & A.H. Jaffe 



(0.8 < n s < 1) than "scale invariance" (n s = 1) gives, a predictable nonzero 
gravity wave contribution with tilt similar to the scalar one, and tiny mean cur- 
vature {yttot ~ 1); (2) same as (1), but with a tiny gravity wave contribution; 
(3) same as (1) but with a subdominant isocurvature component of nearly scale 
invariant tilt (the case in which isocurvature dominates is ruled out); (4) radically 
broken scale invariance with weak to moderate features (ramps, mountains, val- 
leys) in the fluctuation spectrum (strong ones are largely ruled out); (5) radical 
breaking with non-Gaussian features as well; (6) "open" inflation, with quan- 
tum tunneling producing a negatively-curved (hyperbolic) space which inflates, 
but not so much as to flatten the mean curvature (d c ~ (Ha)" 1 , not S> (Ha) -1 , 
where d c = i^ 1 ^^ -1 / 2 ); (7) quantum creation of compact hyperbolic space from 
"nothing" with volume d T which inflates, with dx ~ (Ha) -1 , not S> (Ha) -1 , and 
dx of order d c ; (8) flat (d c = oo) inflating models which are small tori of scale dx 
with dx a few (Ha) -1 in size. It is quite debatable which of the cases beyond (2) 
are more or less plausible, with some claims that (4) is supersymmetry-inspired, 
others that (6) is not as improbable as it sounds. Of course, how likely a priori 
the cases (7) and (8) of most concern to us here is completely unknown, but it 
is the theorists' job to push out the boundaries of the inflation idea and use the 
data to select what is allowed. 

(/) LSS Constraints on the Power Spectrum 

We have always combined CMB and LSS data in our quest for viable models. 
Fig. |l| shows how the two are connected. DMR normalization precisely determines 
<7g for each model considered; comparing with the erg ~ 0.6£7~ T ? 56 target value 
derived from cluster abundance observations severely constrains the cosmological 
parameters defining the models. In Fig. [I], this means the COBE- normalized 
V$(k) must thread the "eye of the needle" in the cluster-band. 

Similar constrictions arise from galaxy-galaxy and cluster-cluster clustering ob- 
servations: the shape of the linear V<s> must match the shape reconstructed from 
the data. The reconstruction shown is from Peacock (1997). The clustering obser- 
vations are roughly compatible with an allowed range 0.15 < T+u s /2 < 0.3, where 
r PS O nr h[O er ./(1.68fi 7 )]- 1 / 2 e-(^(i+^r(2h) 1 /2 ) _ .06) characterizes the density 
transfer function shape. The sCDM model has V ~ 0.5. To get V + zAj/2 in 
the observed range one can: lower h; lower Q nr (ACDM, oCDM); raise 
density parameter in relativistic particles (1.680 7 with 3 species of massless neu- 
trinos and the photons), e.g., as in rCDM, with a decaying v of lifetime and 
r « 1.08O nr h(l + 0.96(m„T d /keV yr ) 2 /3)-i/2. raise q b - tilt Vg < (tCDM), for 
standard CDM parameters, e.g., 0.3 < n s < 0.7 would be required. Adding a hot 
dark matter component gives a power spectrum characterized by more than just 
r. In the post-COBE era, all of these models that lower V + v s /2 have been under 
intense investigation to see which, if any, survive as the data improve. 

(g) Cosmological Radiative Transport 

Cosmological radiative transfer is on a firm theoretical footing. Together with a 
gravity theory (invariably Einstein's general relativity, but the CMB will eventu- 
ally be used as a test of the gravity theory) and the transport theory for the other 
fields and particles present (baryons, hot, warm and cold dark matter, coherent 
fields, i.e., "dynamical" cosmological "constants", etc.), we propagate initial fluc- 
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tuations from the early universe through photon decoupling into the (very) weakly 
nonlinear phase, and predict primary anisotropies, those calculated using either 
linear perturbation theory (e.g., for inflation-generated fluctuations), or, in the 
case of defects, linear response theory. The sources driving their development are 
all proportional to the gravitational potential <3?: the "naive" Sachs- Wolfe effect, 
<I>/3; photon bunching and rarefaction (acoustic oscillations), \^- L , responsible 

for the adiabatic effect and the isocurvature effect; linear-order Thompson 

scattering (Doppler), crrn e v e -g, with ot the Thomson cross section, v e and n e the 
electron velocity and density, and q the photon direction; the (line-of-sight) in- 
tegrated Sachs- Wolfe effect, ~ 2/ Zos there are also subdominant anisotropic 
stress and polarization terms. For primary tensor anisotropies, the sources are 
the two polarization states of gravity waves, hh+ t x] again there are subdominant 
polarization terms. 

Spurred on by the promise of percent-level precision in cosmic parameters from 
CMB satellites (§ [?]), a considerable fraction of the CMB theoretical community 
with Boltzmann transport codes compared their approaches and validated the re- 
sults to ensure percent-level accuracy up to £ ~ 3000 (COMBA 1995). An impor- 
tant goal for COMBA was speed, since the parameter space we wish to constrain 
has many dimensions. Most groups have solved cosmological radiative transport 
by evolving a hierarchy of coupled moment equations, one for each £. Although 
the equations and techniques were in place prior to the COBE discovery for scalar 
modes, and shortly after for tensor modes, to get the high accuracy with speed 
has been somewhat of a challenge. There are alternatives to the moment hierarchy 
for the transport of photons and neutrinos. In particular the entire problem of 
photon transport reduces to integral equations in which the multipoles with I > 2 
are expressed as history-integrals of metric variables, photon-bunching, Doppler 
and polarization sources. The fastest COMB A- validated code, "CMBfast", uses 
this method (Seljak & Zaldarriaga 1996), is publicly available and widely used 
(e.g., to generate some of the power spectra in Fig. ||). 

(h) Secondary Anisotropies 

Although hydrodynamic and radiative processes are expected to play important 
roles around collapsed objects and may bias the galaxy distribution relative to 
the mass (g astrophysics regime in Fig. |l|), a global role in obscuring the early 
universe fluctuations by late time generation on large scales now seems unlikely. 
Not too long ago it seemed perfectly reasonable based on extrapolation from the 
physics of the interstellar medium to the pregalactic and intergalactic medium 
to suppose hydro dynamical amplification of seed cosmic structure could create 
the observed Universe. The strong limits on Compton cooling from the COBE 
FIRAS experiment (Fixsen et al. 1997), in energy 5E Compton coo i/E cmb = Ay < 
6.0 x 10~ 5 (95% CL), constrain the product f exp R 2 exp of filling factor f exp and 
bubble formation scale R exp , to values too small for a purely hydrodynamic origin. 
If supernovae were responsible for the blasts, the accompanying presupernova 
light radiated would have been much in excess of the explosive energy (more 
than a hundred- fold) , leading to much stronger restrictions (e.g., Bond 1996). 

Nonetheless significant "secondary anisotropies" are expected. These include: 
linear weak lensing, dependent on the 2D tidal tensor, a projection of the 3D 

Phil. Trans. R. Soc. Lond. A (1997) 



10 



J.R. Bond & A.H. Jaffe 



tidal tensor d 2 Q /dx l dx J ; the Rees-Sciama effect, 2 f t s <&nl, dependent upon 
the gravitational potential changes associated with nonlinear structure formation; 
nonlinear Thompson scattering, or5n e v e • q, dependent upon the fluctuation in 
the electron density 5n e as well as v e , and responsible for the quadratic-order 
(Vishniac) effect and the "kinematic" Sunyaev-Zeldovich (SZ) effect (moving 
cluster /galaxy effect); the thermal SZ effect, associated with Compton cooling, 
Ii o s , 'pK(x)5(n e T e ), where iPk(x) is a function of x = Ey/T y passing from —2 on 
the Rayleigh- Jeans end to x on the Wein end, with a null at x = 2.83 (i.e., 1863/um 
or 161 GHz); pregalactic or galactic dust emission, ~ f Q s ipdust(xd)Pd, depen- 
dent upon the distribution of the dust density pd and temperature Td through a 
function ifj du st of x d = E 1 /T d . 

Secondary anisotropies may be considered as a nuisance foreground to be sub- 
tracted to get at the primary ones, but they are also invaluable probes of shorter- 
distance aspects of structure formation theories, full of important cosmological 
information. The A:-space range they probe is shown in Fig. |l|. The effect of lens- 
ing is to smooth slightly the Doppler peaks and troughs of Fig. |3|. Q's from 
quadratic nonlinearities in the gas at high redshift are concentrated at high £, 
but for most viable models are expected to be a small contaminant. Thomson 
scattering from gas in moving clusters also has a small effect on Q (although it 
should be measurable in individual clusters). Power spectra for the thermal SZ 
effect from clusters are larger (Bond & Myers 1996); the example in the bottom 
panel of Fig. |3| is for an untilted Hq = 70 COBE-normalized ACDM model, with 
y ~ 2 x 10~ ff (O B h 2 /0.025), still small c.f the FIRAS constraint. (Here and in 
the following, when Hq values are given, the units km s _1 Mpc -1 are implicit.) 

Although C\ may be small, because the power for such non-Gaussian sources 
is concentrated in hot or cold spots the signal is detectable, in fact has been for 
two dozen clusters now at the > 5 sigma level, and indeed the SZ effect will 
soon be usable for cluster-finding. Cg for a typical dusty primeval galaxy model 
is concentrated at higher I associated with galaxy sizes, although a small con- 
tribution associated with clustering extends into the lower I range. These dusty 
anisotropies are now observable with instrumentation on submm telescopes (e.g., 
SCUBA on the James Clerk Maxwell Telescope on Mauna Kea, with the fc-space 
filter shown in Fig. |l|). 



2. CMB Parameter Estimation, Current and Future 

(a) Comparing and Combining CMB Experiments 

We have progressed from the tens of pixels of early AT/T experiments through 
thousands for DMR (Bennett et al. 1996a) and SK95 (Netterfield et al. 1996), 
soon tens of thousands for long duration balloon experiments (LDBs) and even- 
tually millions for the MAP (Bennett et al. 1996b) and Planck (Bersanelli et al. 
1996) satellites. Finding nearly optimal strategies for data projection, compres- 
sion and analysis which will allow us to disentangle the primary anisotropies from 
the Galactic and extragalactic foregrounds and from the secondary anisotropies 
induced by nonlinear effects will be the key to realizing the theoretically-possible 
precision on cosmic parameters and so to determine the winners and losers in 
theory space. Particularly powerful is to combine results from different CMB ex- 
periments and combine these with LSS and other observations. Application of the 
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same techniques to demonstrate self-consistency and cross-consistency of experi- 
mental results is essential for validating conclusions drawn from the end-product 
of data analysis, e.g., the power spectra in bands as shown in Fig. || and the 
cosmic parameters they imply. 

Current band-powers are shown in the upper panel of Fig. [2|. The first lesson 
of Figs. [2|,|| is that, in broad brush stroke, smaller angle CMB data (e.g., SP94, 
SK95, MSAM, MAX) are consistent with COBE- normalized C^'s for the untilted 
inflation-based models. It is possible that some of the results may still include 
residual contamination, but it is encouraging that completely different experi- 
ments with differing frequency coverage are highly correlated and give similar 
bandpowers, e.g. DMR and FIRS (e.g., Bond 1996), SK95 and MSAM (Netter- 
field et al. 1996, Knox et al. 1998). Lower panels compress the information into 
9 optimal bandpower estimates derived from all of the current data (see Bond, 
Jaffe & Knox 1998b for techniques). 

The few data points below i < 20 are mainly from COBE's DMR experiment. 
Clearly the Grange spanned by DMR is not large enough to fix well the cos- 
mological parameter variations shown in the right panels, but combining CMB 
anisotropy experiments probing different ranges in £-space improves parameter 
estimates enormously because of the much extended baseline: it is evident that 
n s can be reasonably well determined, low £1 open models violate the data, but 
Q\ cannot be well determined by the CMB alone. 

(b) DMR and Constraints on Ultra-large Scale Structure 

DMR is fundamental to analyses of the VLSS region and ULSS region, and is 
the data set that is the most robust at the current time. The average noise in 
the 53+90+31 GHz map is about 20/i-fT per fwhm beam (~ 7°), and there are 
about 700 of these resolution elements outside of the Galactic disk cut (about 
4000 2.6° DMR pixels with 60 fiK noise). The signal is about 37 \xK per beam: 
i.e., there is a healthy signal-to-noise ratio. The signal-to-noise for widespread 
modes (e.g., multipoles with I < 15) is even better. Indeed, even with the much 
higher precision MAP and Planck experiments we do not expect to improve the 
results on the COBE angular scales greatly because the 4-year COBE data has 
sufficiently large signal-to-noise that one is almost in the cosmic variance error 
limit (due to realization to realization fluctuations of skies in the theories) which 
cannot be improved upon no matter how low the experimental noise. 

Wiener-filtered maps shown in Fig. ||] give the statistically-averaged signal given 
the data and a best-fit signal model. These optimally-filtered maps are insensi- 
tive to modest variations in the assumed theory. The robustness of features in 
the maps as a function of frequency and the weak frequency dependence in the 
bandpowers are strong arguments that what is observed is on the sky with a pri- 
mary anisotropy origin, made stronger by the compatible amplitudes and positive 
cross-correlations with the FIRS and Tenerife data sets. 

Recall the "beyond our horizon" land in Fig. [l] is actually partly accessible 
because long waves contribute gentle gradients to our observables. The DMR data 
is well suited to probe this regime. Constraints on such "global parameters" as 
average curvature from COBE are not very good. Obviously it is much preferred 
to use the smaller angle data on the acoustic peak positions. The COBE data can 
be used to test whether radical broken scale invariance results in a huge excess or 
deficit of power in the COBE /c-space band, e.g., just beyond k~ 1 ~ Hq 1 , but this 
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has not been much explored. The remarkable non- Gaussian structure predicted 
by stochastic inflation theory would likely be too far beyond our horizon for its 
influence to be felt. The bubble boundary in hyperbolic inflation models may be 
closer and its influence shortly after quantum tunneling occurred could possibly 
have observable consequences for the CMB. Theorists have also constrained the 
scale of topology in simple models (Fig. |4|). Bond, Pogosyan &; Souradeep (1997, 
1998) find the torus scale is d T /2 > 1.1{2Hq 1 ) = 6600 b" 1 Mpc from DMR 
for flat equal-sided 3-tori at the 95% confidence limit, slightly better than other 
groups find since full map statistics were used. The constraint is not as strong 
if the repetition directions are asymmetric, > 0.7(2i7 _1 ) for 1-tori from DMR. 
It is also not as strong if more general topologies are considered, e.g., the large 
class of compact hyperbolic topologies (Bond, Pogosyan & Souradeep 1997, 1998, 
Cornish et al. 1996, 1998, Levin et al. 1997, 1998). 

(c) Cosmic Parameters from All Current CMB Data 

We have undertaken full Bayesian statistical analysis of the 4 year DMR (Ben- 
net et al. 1996a), SK94-95 (Netterfield et al. 1996) and SP94 (Gundersen et al. 
1995) data sets, taking into account all correlations among pixels in the data and 
theory (Bond & Jaffe 1997). Other experiments available up to March 1998 were 
included by using their bandpowers as independent points with the Gaussian er- 
rors shown in Fig. |2[ We have shown this approximate method works reasonably 
well by comparing results derived for DMR+SP94+SK95 with the full analysis 
with those using just their bandpowers (Jaffe, Knox & Bond 1997). We have also 
shown that a significant improvement in accuracy is possible if instead of the 

1/2 

average and 1 sigma limits on the experimental bandpowers or on Cg , one 
uses ln(Ce + xb), where xb is related to the noise of the experiment (Bond, Jaffe 
& Knox 1998b). This includes some of the major non-Gaussian deviations in the 
bandpower likelihood functions; results using this more accurate approach, and 
incorporating the very recent CAT98 and QMAP data, will be reported elsewhere 
(Bond, Jaffe &: Knox 1998c, in preparation). Other groups have also calculated 

1/2 

parameter constraints using the C B bandpower approach (e.g., Lineweaver &: 
Barbosa 1997, Hancock & Rocha 1997, Lineweaver 1998). 

With current errors on the data, simultaneously exploring the entire parameter 
space of § [T] is not useful, so we restricted our attention to various subregions of 
{n B h 2 , O cdm h 2 , Mm h 2 , Q fc h 2 , ^Ah 2 , v s , v t , a s }, such as {a 8 , n s , h | fixed t , rj s h 2 }, 
where 0^=0 and 17a is a function of hio or f2A=0 and fife > is a function of 
hio- The age of the Universe, to, was chosen to be 11, 13 or 15 Gyrs. A recent 
estimate for globular cluster ages with the Hipparcos correction is 11.5 ± 1.3 Gyr 
(Chaboyer et al. 1998), with perhaps another Gyr to be added associated with 
the delay in globular cluster formation, so 13 Gyr is a good example. We consid- 
ered the ranges 0.5 < n s < 1.5, 0.43 < h < 1, and 0.003 < B h 2 < 0.05. The old 
"standard" nucleosynthesis estimate was n B h 2 = 0.0125, but the preferred one 
is now 0.025. We assumed reheating occurred sufficiently late to have a negligible 
effect on Cg, although this is by no means clear. C/s for sample restricted param- 
eter sequences are shown in Fig. We made use of signal-to-noise compression of 
the data (by factors of 3) in order to make the calculations of likelihood functions 
such as C(as, n s , h | fixed to, O^h 2 ) more tractable (without loss of information or 
accuracy) . 
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Figure 4. The first column shows unfiltered 140° diameter dmr A+B maps centered on the 
North Galactic Pole, the second shows them after Wiener-filtering (with monopole, dipole and 
quadrupole removed), the third the South Pole version, with the nth contour as noted and 
negative contours heavier than positive ones. The Wiener maps use a model which fits the 
correlation function and amplitude of the DMR data (specifically, the n 3 — l sCDM model was 
used, but insensitive to even rather significant variations). The maps have been smoothed by a 
1.7° Gaussian filter, all is 53+90+31t4+_B. Although higher noise results in filtering on greater 
angular scales, the large scale features of all maps are the same. This is also borne out by detailed 
statistical comparisons map to map. The last column shows some theoretical realizations, after 
optimal filtering. The first two rows are the NGP and SGP for a n 3 =l CDM model. The lower 
two rows are for a 3-torus topology, with repetition length dr = 9000 h _1 Mpc, 1.5 times the 
horizon radius, in all three directions, a model strongly ruled out because of the high degree of 
positive correlation between the North and Southern hemispheres that the periodicity induces 
(Bond, Pogosyan & Souradeep f998). Highly correlated patterns also exist for small compact 
hyperbolic models and lead to constraints on manifold size. 

The n s constraints are quite good. If a$ is marginalized for the tilted ACDM 
sequence with Hq=50, with DMR only the primordial index is n s = 1.02^25 

with no gravity waves and fi=0, and 1.02^g with gravity waves and vt=v s , 
rather encouraging for the nearly scale invariant models preferred by inflation 
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theory. Because the gravitational potential changes at late time with A^O, the 
integrated Sachs- Wolfe effect gives more power in Ci at small £, so the preferred 
n s steepens to compensate. When A is marginalized in the 13 Gyr tilted ACDM 
sequence, n s = 1.17 ±0.31 is obtained. For this sequence, when all of the current 
CMB data are used we get L02j/gf for Hq = 50 (and tt A = 0, the tilted sCDM 
model sequence) and l.OOt'^ for Hq = 70 (and = 0.66). Marginalizing over 
Hq (i.e., A) gives 1.0lj? 04 with gravity waves included, 0.98l; g if they are not. 

The marginalized 13 Gyr tilted oCDM sequence gives l.OOl'ojj. 

Hq and Oa for fixed age are not that well determined by the CMB data alone, 
as can be seen from the dotted lines in Fig. After marginalizing over all n s , we 
get Hq < 75 at la, but effectively no constraint at 2a. The strong dependence of 
the position of the acoustic peaks on means that the oCDM sequence is better 
restricted: VL to t ~ -7 is preferred; for the 13 Gyr sequence this gives Hq sa 53 and 
for the 11 Gyr sequence Hq m 65. 

Calculations of defect models (e.g. strings and textures) give Q's that do not 
have the prominent peak that the data seem to indicate (Pen, Seljak & Turok 
1997, Allen et al. 1997). 

( d ) Cosmic Parameters from Current LSS plus CMB Data 

Combining LSS and CMB data gives more powerful discrimination among the 
theories, as Fig. |l] illustrates visually and Fig. |^ shows quantitatively. The ap- 
proach we use here and in Bond & Jaffe (1997) to add LSS information to the 
CMB likelihood functions is to design prior probabilities for T + u s /2 and crgf^',? 6 , 
reflecting the current observations, but with flexible and generous non-Gaussian 
and asymmetric forms to ensure the priors can encompass possible systematic 
problems in the LSS data. For example, our choice for ogfi^;? 6 was relatively flat 
over the 0.45 to 0.65 range. (Explicitly we used 0.55j^oft oi' with the two error 
bars giving a Gaussian and a top hat error so that the net result is generously 
flat over the total ±la range. For T + u s /2, we used Using the Pea- 

cock 1997 reconstructed linear power spectrum shown in Fig. |l] would give more 
stringent constraints for the shape e.g., Gawiser & Silk 1998.) 

Using all of the current CMB data and the LSS priors, for the 13 Gyr ACDM 
sequence with gravity waves included, we get n s = 1.00l;°| and Hq = 72 ± 3 
(Oa ~ 0.7), respectively, when Hq and n s are marginalized; with no gravity waves, 
0.96±;g| and Hq = 70 ±3 are obtained; and for an AhCDM sequence, with a fixed 
ratio ^Ihdm/^nr = 0.2 for two degenerate massive neutrino species, n s rj 0.97^;q2 
and Hq ps 57tl are obtained, revealing a slight preference for £7a ~ 0.3. For the 
15 Gyr ACDM sequence, the tilts remain nearly scale invariant and near 0.6: 
0.98±;{$ and Hq = 57 ± 3 (Q A « 0.6) with gravity waves, 0.95± g| and 54 ± 3 
(Q\ ss 0.5) without. 

For the 13 Gyr oCDM sequence, the likelihood peak for the CMB+LSS data 
is shifted relative to using the CMB data alone because the best fit CMB-only 
models have as too low compared with the cluster abundance requirements. Al- 
though the Hq sa 54+J value (fltot ~ 0.6) is close to the CMB-only one, the 
maximum likelihood is significantly below the ACDM one. Hq is larger for the 11 
Gyr oCDM sequence, but £ltot is about the same, and the likelihood is still low. 

Should these small error bars be taken seriously? It seems unlikely that as from 
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Figure 5. Likelihood curves for fixed-age ACDM and oCDM sequences, marginalized over as 
and n s . A Gaussian approximation to the likelihood places 1,2,3a at the horizontal dashed 
lines. The la ranges are explicitly given in the text. The dotted curves are for CMB only, solid 
for CMB+LSS. The right panels are equivalent to the left, but translated to Q nr (1 — Qa for 
ACDM, fitot for oCDM). The curves shown are for no GW, but there is little difference if GW are 
included. The sequence labelled m v is for the 13 Gyr AhCDM sequence with a fixed !J m i//n„ r 
ratio of 0.2, and 2 degenerate neutrino species. Even this case favours slightly a nonzero Qa- 
In the upper left panels, the 68% confidence limits for the values of Q. m estimated using the 
luminosity-distance relation for Type I supernovae from Perlmutter et al. and Craig et al. are 
also shown. Note that the CMB data alone slightly favours a value of Q.tot < !■• The absolute 
likelihood for the CMB+LSS data strongly favours the ACDM over the oCDM sequences. When 
just the fully analyzed DMR+SK95+SP94 data are used with LSS, the marginalized results 
for ACDM are remarkably similar, with the long baseline between DMR and SK95 fixing the 
freedom in n s (although values of n 3 about 1.1 are now preferred). When only DMR is used 
along with LSS, n s is not nailed down, and the resulting freedom implies f2„ r = 1 models are 
not disfavoured. The horizontal error bars in the upper right panel show the la range of finr 
for Q nr + Qa = 1 models inferred from the supernova la observations of Perlmutter et al. 1998 
(upper) and Reiss et al. 1998 (lower). 
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cluster abundances will change much; and, as we have seen, the DMR results are 
quite robust. Although largely driven by just the DMR plus LSS results, the 
smaller angle CMB results lock in the tilt, and as the CMB data improves some 
adjustment might occur, but not a drastic one unless we have made a major 
misinterpretation in the nature of the CMB signals observed at intermediate 
angles. 

If we were to marginalize over to as well, it is clear that Hq would not be as 
well determined, but n s and either fij. or il\ would be. If the parameter space 
is made even larger, near degeneracies among some cosmic parameters become 
important for CMB data alone, and these are only partially lifted by the LSS data 
(e.g., Efstathiou & Bond 1998). In particular, this restricts the ultimate accuracy 
that can be achieved in the simultaneous determination of and This will 
become an issue when the quality of the CMB data improves, as described in the 
next subsection, but for now one must bear in mind the constrained space used 
when interpreting the current precision quoted on parameter estimation. 

( e ) Cosmic Parameters from the CMB Future 

The expected error bars on the power spectrum from MAP and Planck (Bond, 
Jaffe & Knox 1998a, Bond, Efstathiou & Tegmark 1997, hereafter BET) shown in 
Fig. |2] illustrate that even quite small differences in the theoretical C/s and thus 
the parameters can be distinguished. Quite an industry has developed forecasting 
how well future balloon experiments (Maxima, Boomerang, ACE, Beast, Top 
Hat), interferometers (VSA, CBI, VCA) and especially the satellites MAP and 
Planck could do in measuring the radiation power spectrum and cosmological 
parameters if foreground contamination is ignored (Knox 1995, Jungman et al. 
1996, BET, Zaldarriaga, Spergel, & Seljak 1997, White, Carlstrom & Dragovan 
1997). Forecasts like these were quite influential in making the case for MAP and 
Planck. 

Table || gives some examples of what can be obtained using only CMB data 
(BET). The experimental parameters chosen are given in Table [[]. The durations 
chosen were appropriate for the types of experiments, e.g., about a week for long 
duration balloon experiments and about two years for satellite experiments. The 
temperature anisotropies were assumed to be Gaussian-distributed, and among 
the > 17 parameters of § En a restricted 9 parameter space was used: 5 densities, 
Q nr , Qhdm, ^fcj ^A}h , the Compton depth tc, the scalar tilt, n s , the total 
bandpower for the experiment (Q) B in place of V$(k n ), and the ratio of tensor 

(T) ( S) 

to scalar quadrupole powers, r ts = C\ jC 2 , in place of v t . Just like Vgw/'P^, 
r ts is a sensitive function of vt, but also depends on v s — u t , 0,a, etc. (Bond 1996). 
In this space, recall that h 2 = £.,(%h 2 ) is a dependent quantity. 

Except for the integrated Sachs- Wolfe effect at low £, the angular pattern of 
CMB anisotropies now is a direct map of the projected spatial pattern at red- 
shift ~ 100, dependent upon the cosmological angle-distance relation, which is 
constant along a line relating 0/%h 2 and Q\h 2 for fixed Q nr h 2 . This defines a 
near-degeneracy between f2/% and Q\ broken only at low t where the large cos- 
mic variance precludes accurate determination of both parameters simultaneously 
(e.g. BET, Zaldarriaga, Spergel, & Seljak 1997, Efstathiou & Bond 1998, Eisen- 
stein, Hu & Tegmark 1998). Other cosmological observables are needed to break 
this degeneracy. A good example is Type I supernovae. If they are assumed to 
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Table 1. The noise is per fwhm pixel, in fiK. Max is Maxima, TH is TopHat, and Bst is Beast, 
which also has a 30 GHz channel (26', 53/^K) included in the analysis. Bm is Boomerang, which 
also has 220 and 430 GHz channels not included in the analysis. The North America test flights 
covered .008 and .003 of the sky, and had fwhm of 35' and 19', respectively. Boost is a bolometer 
based version of Beast that is an example of an ultra-long duration balloon experiment, of order 
100 days, with channels at 100 GHz (9', 17/xK), 140 GHz(6', 24/iK), (4', 37^K), covering .07 of 
the sky, with l cu t ~ 6. The Planck numbers are for the HFI channels; a 350 GHz (4', 43/iK) was 
also included in the analysis. See BET for the 3 channels of LFI used. 
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Table 2. Sample idealized MAP and Planck parameter error forecasts, for a 9 parameter inflation 
family of models, with standard CDM the "Target Model". See BET for methods. ^Ah 2 is deter- 
mined with S7feh 2 fixed, and fifch 2 is determined with fiAh 2 fixed, because of the angle-distance 
near-degeneracy (e.g., Efstathiou & Bond 1998); the other parameters are insensitive to fixing 
either, or neither. The ranges for Ho, fish 2 are absolute, but the errors are relative ones. The 
forecasted errors obviously represent a great leap forward from current errors and from what 
is conceivable with non-CMB probes. Amplitude parameters are highly correlated with rc, but 
this can be partly broken when other information is included, e.g., on the abundance of clus- 
ters. The third column is an optimistic forecast of what one can do with balloons by combining 
MAXIMA, TopHat, Boomerang and BeastI with DMR (see Fig. ||). TopHat, Boomerang, and 
BeastI would be long duration balloon flights, lasting about a week over the Antarctic. The 
parameters used are given in Table |e|. It is unclear that systematics will be sufficiently small for 
the LDB experiments to fulfill this promise. 
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be "standard candles", then their degeneracy is along lines of equal luminosity- 
distance, which is sufficiently different from the equal angle-distance lines to allow 
good separate determination. 

If the polarization power spectrum can be measured with reasonable accuracy, 
errors on some parameter such as rj s would improve (Zaldarriaga et al. 1997). 
However the polarization power spectrum is about a hundred times lower than 
the total anisotropy, and the gravity wave induced polarization is substantially 
tinier than this at the low £ needed for rt s improvement. We do not know if the 
foreground polarization will hopelessly swamp this signal. 

Error forecasts do depend upon the correct underlying theory. In Table ||, un- 
tilted sCDM was chosen as the target model, but the values shown are indicative 
of what is obtained with other targets (BET). The third column gives errors fore- 
casted for balloon experiments, the bolometer-based TopHat, Boomerang, and 
MAXIMA and the HEMT-based BEAST. (URLs to home pages are given in the 
references.) ^-cuts were included to reflect the limited sky coverage these experi- 
ments will have. Adding DMR to extend the ^-baseline diminishes the forecasted 
errors. 

We adopt the current beam sizes and sensitivities for MAP and Planck used in 
BET, improvements over the original proposal values. Of the 5 HEMT channels 
for MAP, BET assumed the 3 highest frequency channels, at 40, 60 and 90 GHz, 
will be dominated by the primary cosmological signal (with 30 and 22 GHz chan- 
nels partly contaminated by bremsstrahlung and synchrotron emission). MAP 
also assumes 2 years of observing. For Planck, BET used 14 months of observing, 
the 100, 65, 44 GHz channels for the HEMT-based LFI (but not the 30 GHz chan- 
nel), and the 100, 150, 220 and 350 GHz channels for the bolometer-based HFI 
(but not the dust-monitoring 550 and 850 GHz channels). The highest resolution 
for MAP is 13' fwhm, the highest for Planck is 4'. 

These idealized error forecasts do not take into account the cost of separating 
the many components expected in the data, in particular Galactic and extragalac- 
tic foregrounds, but there is currently optimism that the Galactic foregrounds at 
least may not be a severe problem e.g., (Bersanelli et al. 1996), although low fre- 
quency emission near 100 GHz by small spinning dust grains (Leitch et al. 1997, 
Draine &; Lazarian 1998) may emerge as a new significant source. There is more 
uncertainty about the extragalactic contributions in the submm and radio. 

Although we may forecast wonderfully precise power spectra and cosmic pa- 
rameters for the simplest inflation models in Table once we consider the more 
baroque models with multifeatured spectra the precision drops(e.g., Souradeep et 
al. 1998). Given that all of our CMB and LSS observations actually access only 
a very small region of the inflation potential, imposing theoretical "prior" costs 
on highly exotic post-inflation shapes over the observable bands is reasonable. 
Nonetheless, if the phenomenology ultimately does teach us that non-baroque in- 
flation and defect models fail, the CMB and LSS data will be essential for guiding 
us to a new theory of fluctuation generation. 

We would like to thank George Efstathiou, Lloyd Knox, Dmitry Pogosyan, 
and Tarun Souradeep for enjoyable collaborations on a number of the projects 
highlighted in the text. 
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